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We study a dense colloidal suspension confined between two quasiparallel glass plates as a model 
system for a supercooled liquid in confined geometries. We directly observe the three-dimensional 
Brownian motion of the colloidal particles using laser scanning confocal microscopy. The particles 
form dense layers along the walls, but crystallization is avoided as we use a mixture of two particle 
sizes. A normally liquid-like sample, when confined, exhibits slower diffusive motion. Particle 
rearrangements are spatially heterogeneous, and the shapes of the rearranging regions are strongly 
influenced by the layering. These rearranging regions become more planar upon confinement. The 
wall-induced layers and changing character of the spatially heterogeneous dynamics appear strongly 
connected to the confinement induced glassiness. 

PACS numbers: 64.70.pv, 61.43.Fs, 82.70.Dd 



I. INTRODUCTION 

As a glass-forming liquid is cooled, its viscosity in- 
creases smoothly but dramatically by many orders of 
magnitude. The macroscopic divergence in viscosity is 
related to the divergence in the microscopic structural re- 
laxation time, or a-relaxation time. A conceptual expla- 
nation is the Adams and Gibbs hypothesis, which states 
that the flow in a supercooled liquid involves the cooper- 
ative motion of molecules and that the structural arrest 
at the glass transition is due to a divergence of the size 
of these cooperatively rearranging regions (CRRs) [1]. 

Computer simulations and experiments have explored 
the sizes and shapes of regions of cooperatively mov- 
ing molecules as a liquid's glass transition is approached 
[2, 3]. A direct means of probing the dynamic length 
scales of glass-forming liquids is by confining them to 
smaller volumes, such as within thin films and nanopores. 
Confinement can either increase, decrease or even main- 
tain a material's glass transition temperature Tq [4, 5]. 
Both simulation and experiment suggest that the effect 
on Tq depends on the nature of the interaction between 
the sample and its confining boundary [6-12]. Attrac- 
tive interactions may result in an increase in Tq whereas 
repulsive interactions may result in a decrease [10, 13]. 
Frustration of structural ordering, via a rough surface for 
example, can also play a key role, although this can ei- 
ther cause slower or faster dynamics [7, 11, 14]. Whether 
or not the restriction of the length scales accessible to 
CRRs is responsible for the variation in Tq remains to 
be seen due to the inability to directly observe molecular 
interactions within glass-forming liquids. 

Instead of studying molecular glass-formers, we 
use dense colloidal suspensions of sterically-stabilized 
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micrometer-sized spherical particles. Colloidal suspen- 
sions have often been used as experimental models of 
a hard sphere glass [15, 16]. We confine our samples 
within a planar volume formed by two quasi-parallel 
solid surfaces [17], similar to confined colloids studied 
by other groups [11, 18, 19]. We use high-speed con- 
focal microscopy to rapidly visualize and acquire three- 
dimensional images of the particle positions [20, 21] . Sub- 
sequent image analysis lets us track the individual par- 
ticle trajectories, providing an accurate picture of the 
cooperatively rearranging groups of particles. Near the 
colloidal glass transition (0g « 0.58 [15]), particles rear- 
range in groups characterized by a length scale of ~3-6 
particle diameters [20, 22]. 

In this manuscript we further investigate our results 
from prior experiments that studied confined samples, 
as pictured in Fig. 1 [17]. Here we focus specifically 
on the nature of cooperative rearrangements within the 
confined sample and how they relate to the system's in- 
creased glassiness. In these experiments we found that 
confinement induces glassy behavior at concentrations in 
which the bulk behavior is still liquid-like. Here, we show 
that confining colloidal liquids within this planar volume 
results in cooperatively rearranging groups of particles 
that are similarly planar shaped. We find that the flat- 
tening shapes of the cooperatively rearranging groups are 
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FIG. 1: (Color online) Sketch of sample chamber (not to 
scale). The small particles are 1.18 ^m in radius and are 
shaded to indicate their fluorescent dye. The large particles 
are 1.55 fim in radius and drawn in white to indicate their 
lack of dye, making them invisible to the confocal microscope. 
One of the boundaries is a coverslip, rather than a glass slide, 
indicated by the thinner line. 
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correlated with the overall slowing of the dynamics, sug- 
gesting a connection between confinement, wall- induced 
structure, and glassy behavior. 

Understanding the effects of confinement on the glass 
transition may help us understand the glass transition in 
the bulk. Perhaps more importantly, understanding the 
properties of confined fluids also has direct relevance with 
lubrication [23], the flow of liquids through microfluidic 
devices [24, 25], and the kinetics of protein folding [26]. 

II. EXPERIMENTAL DETAILS 
A. Colloidal samples and microscopy 

We use spherical colloidal poly-methyl-methacrylate 
(PMMA) particles that are sterically stabilized to pre- 
vent interparticle attraction [15, 21]. The particles are 
suspended in a mixture of solvents, cyclohexylbromide 
and cis- and trans-decalin, which matches both their den- 
sity and index of refraction [21]. While our sample is 
similar to other types of colloidal suspensions that act 
like hard spheres [15], the cyclohexylbromide in our sol- 
vent mixture induces a slight charge on the surfaces of 
the particles. Thus, the particles have a slightly soft re- 
pulsive interaction in addition to their hard sphere core. 
To prevent crystallization, which would be readily in- 
duced by the smooth walls in our thin planar geometry 
[27, 28], we use a binary mixture of particles with hydro- 
dynamic radii of asmaii = 1-18 /im and aiargc = 1-55 /im. 
The number ratio is approximately Ng/N]^ = 3.5, and 
the individual volume fractions are approximately (ps = 
0.26,01, = 0.16, so the total overall volume fraction is 
(j) = 0.42 ± 0.05. The uncertainty of (j) arises from the 
difficulty in precisely determining the individual species' 
particle size, the polydispersity of particle sizes (~ 5% 
for both species), and difficulties in determining the rel- 
ative volume fractions of the two species [29] . A study of 
a similar colloidal mixture found the glass transition for 
bulk samples to be at 0^ ss 0.58 [30]. 

We use laser scanning confocal microscopy to view the 
sample [21]. We can acquire a three-dimensional image of 
the sample by scanning a 50 x 50 x 20 fim"^ region (equal 
to 256 X 256 x 100 pixels). We use Visitech's "vt-Eye" 
confocal system which can scan this volume in 2.0 sec- 
onds. This is much faster than the time for particles to 
diffuse their own diameter, which is ~100 seconds in our 
samples. We acquire sequences of three-dimensional (3D) 
confocal images every 2.0 seconds for up to 45 minutes. 
By scanning different locations, we observe the behavior 
at different chamber thicknesses ranging from ^6 /im to 
~19 /im in addition to the sample's bulk. Data repre- 
senting the 'bulk' of our sample is acquired from a 20 /im 
thick subvolume in the thicker region of the sample cham- 
ber that is over 15 /im away from the chamber's walls in 
order to avoid any boundary effects. 

The small particles are dyed with Rhodamine dye [21] 
and the larger ones are left undyed. Thus the data in our 




FIG. 2: Typical 2D confocal microscope image showing par- 
ticles immediately adjacent to one of the chamber walls. The 
circled particles are stuck to the glass, and the others move 
freely. There are also undyed particles also stuck to the sur- 
face, as well as undyed mobile particles, which are not visible 
in this confocal image. The scale bar indicates 10 /xm. 

results are for the smaller particles only. Each image is 
post-processed to find particle positions with an accuracy 
of 0.05 /im in x and y (parallel to the walls) and 0.1 fj,m 
in z (perpendicular to the walls, and parallel to the opti- 
cal axis of the microscope). Given that the particles do 
not move much between images, we can link the particle 
positions in time to get 3D trajectories of the particles' 
motion throughout the sample volume [21, 31]. 

B. Sample chambers 

Our goal is to study our sample with a range of confine- 
ment thicknesses. Here we focus on "thin film like" con- 
finement between two fiat surfaces. We achieve this by 
constructing a wedge shaped sample chamber, as shown 
in Fig. 1. We build the chamber using a glass slide, a rect- 
angular glass coverslip, and a narrow piece of a ~60 /xm 
thick Mylar film, employing a method similar to the one 
used by Refs. [32, 33]. Using UV-curing epoxy (Norland 
68) we attach the Mylar film near one end of one side of 
the glass slide so that it runs perpendicular to the slide's 
length. Next, the glass coverslip is laid across the slide so 
that one end is raised up by the Mylar film. Meanwhile, 
the coverslip's opposite end is clamped down, ensuring 
the thinnest gap size possible. We seal the sample cham- 
ber shut with epoxy, except for two small air holes; the 
sample is added via one while air escapes via the other. 
After adding the sample, the two openings are sealed with 
epoxy. The chamber's shape is described in Fig. 1: a very 
long chamber with a broad range of thicknesses. Due to 
the Mylar film, the glass surfaces are not parallel but very 
slightly angled at 0.4° relative to one another. Within our 
field of view, the change in our sample's thickness due to 
our sample chamber's slight taper is less than 0.3 /im, 
which is negligible for all but the thinnest regions. We 
do not see any influence of the taper in any of our results 
(discussed further below), suggesting it is reasonable to 
consider the two boundaries as locally quasi-parallel. We 
define y as the direction along which H varies. 

When we fill our slides with sample, a small fraction of 
particles stick to the sample chamber's walls. Typically 
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less than 20% of the walls' area is coated with stuck par- 
ticles [34]. The stuck particles are easy to identify as 
their apparent motion, due to noise inherent to particle 
tracking, is much less than the other particles. An image 
showing the locations of some stuck particles is shown in 
Fig. 2. Other observations confirm that both large and 
small particles stick to the walls [34]. We find that the 
particles stick to the surfaces of the glass slides only dur- 
ing the initial loading of the sample chamber with colloid. 
The stuck particles remain stuck indefinitely, through a 
van der Waals attraction to the glass, and are a perma- 
nent feature of the surface. The mobile particles do not 
stick to the sample's glass boundaries over time - during 
the experiments they never are seen to stick, and over 
several months the amount of particles stuck to the glass 
does not appear to change. In fact, the mobile parti- 
cles are repelled from the glass boundaries by a relatively 
weak Coulombic interaction; in other words, during the 
course of the experiment, the only particle- wall interac- 
tion is a weakly repulsive one. In a sample of dilute col- 
loids, we observe that the concentration of particles is low 
at the wall and approaches the bulk value quickly, within 
0.5 ±0.1 fim, suggesting that the Debye screening length 
is « 0.4 /j,m at most and more likely « 0.2 — 0.3 /xm, in 
agreement with prior observations [35]. The stuck parti- 
cles are expected to slightly slow adjacent particles [36], 
which has been confirmed in our experimental data [34]. 

Particles do interact with the wall hydrodynamically. 
In the same dilute suspension, we measure particle mobil- 
ity near the glass walls, with measured diffusivity shown 
in Fig. 3 as a function of the distance z from the wall. The 
behavior (symbols) is in good agreement with Faxen's 
Law (dashed line) [37, 38] which quantifies the hydro- 
dynamic influence of a planar boundary. Of course, the 
hydrodynamic behavior is modified in confinement ap- 
proaching quasi-two-dimensional situations, where the 
sample chamber thickness H is comparable to the parti- 
cle size 2a [39-41]. We do not consider experiments that 
are this thin; our observations all have H > 6a. More sig- 
nificantly, for the larger volume fractions we consider in 
this work, the hydrodynamic interaction will be screened 
by the other particles, and so will not depend so strongly 
on the distance from the wall [41, 42]. 



III. RESULTS 

A. Wall-induced structure 

We use the positions of the stuck particles to measure 
the local thickness of the sample chamber. To do this we 
find the number density n{z) as a function of the distance 
z between the walls, shown in Fig. 4 for (a) the mobile 
particles and (b) the stuck particles. The maximum of 
each peak in (b) corresponds to the approximate position 
of the centers of the small particles stuck to the sample's 
walls. These positions are marked by the vertical dashed 
lines in Fig. 4, whose separation indicate the effective 



1.50 r 

1.25 
1.00 
0.75 
0.50 
0.25 
0.00 



z (jJsn) 

-2 2 4 6 8 10 12 14 16 



A D 



covers lip 



2 4 6 8 10 12 14 
da 



FIG. 3: (Color online) Measurements of local diffusion con- 
stants as a function of the distance z to the wall, normalized 
by the small particle radius a — 1.18 fim. Inset: sketch in- 
dicating that 2 = corresponds to the particle touching the 
wall. 
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FIG. 4: (Color online) The number density n as a function of 
the distance z between the walls, for the visible (small) parti- 
cles, (a) All mobile particles, (b) All immobile particles. The 
vertical lines in both indicate the position of the centers of the 
visible particles stuck to the walls. For this data, the distance 
between the two positions is H = 6.25 fim, the effective local 
chamber thickness. 
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FIG. 5: (Color online) Number density n as a function of the 
distance z between the walls. The thickness H is as labeled 
(in microns). The values of H in terms of Osmaii are 5.30, 5.62, 
6.83, 6.86, 10.7, and 13.4. The curves are vertically offset for 
clarity, where the offset is proportional to H. Where there is 
an asymmetry in the height of the z; « peak and the z k, H 
peak, it is due to one wall having more stuck particles on it, 
thus decreasing the room available for mobile particles. 



local chamber thickness H. Since only the small particles 
are visible to the microscope, the actual thickness is iJ + 
2asmaii = H + 2.36 fim. The mean particle radii are 
known only to within ±0.02 ^m, while our uncertainty 
in their z positions is 0.1 fim. By averaging over tens of 
stuck particles we can determine H to within 0.01 fim. 

Figure 4(a) shows layering of particles near the sam- 
ple walls, which has been seen in both computer simula- 
tions [6, 43] and experiments [11, 28, 32, 36]. Comparing 
Fig. 4(a) to (b) we see that the boundary layers of the 
mobile particles are offset from those of the stuck parti- 
cles. The offset is due to Coulombic repulsion between 
the glass walls and PMMA particles, and is about 0.4 /im 
in all cases. Using differential interference contrast (DIG) 
microscopy, we confirm that the large particles also form 
layers, albeit in positions shifted due to their size. Our re- 
sults are qualitatively in agreement with simulations that 
studied layering of binary mixtures of particles near walls 
[44, 45], and are fairly similar to observations of layering 
in single-component colloidal samples [11, 28, 36]. 

Figure 5 displays the way layering changes with H. 
The peaks of n{z) are tallest and thinnest next to the 
walls. Subsequent layers are shorter and wider, pre- 
sumably as the correlations between particle positions 
become diluted through the presence of two particle 
sizes [44]. Note that we do not see any "quantiza- 
tion" effects for particular values of H [46]. For ex- 
ample, some packing effects were seen in simulations at 
H = 2masmaii + 2naiargc for integer values to, n, but these 
effects are too subtle to be resolved given the relatively 
few values of H for which we have experimental data [44] . 
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FIG. 6: (Color online) (a) The mean square displacement 
for our sample over a range of thicknesses [17]. The dashed 
line has a slope of 1.0. (b) Plots of the corresponding non- 
Gaussian parameter for each thickness. The x superscript of 
af is to indicate that the non-Gaussian parameter is only cal- 
culated using the x displacements (parallel to the wall, and 
perpendicular to the slight gradient in H). The inset is a mag- 
nification of the curves for H > 15.8 /im, with each curve's 
local maxima labeled, corresponding with At* for the data 
at these thicknesses, (c) Components of the MSD curves. 
Light gray (red) curves are the x- and y-components of mo- 
tion (parallel to the walls) and the dark gray (blue) are the 
z-component of motion (perpendicular). 



B. Sample-averaged dynamics 

Before we consider the specific infiuence of the particle 
layers on the particle motion, we will quantify the average 
motion of the sample. This is done by calculating the 
mean square displacement (MSD) as 



where the average is taken over all particles i and all ini- 
tial times t. Analogous formulas apply for (j/^) and (z^). 
Figure 6(a) shows that the motion parallel to the walls 
slows dramatically with confinement (decreasing H, as 
indicated). For values less than iJ « 16 /xm « 14 agmaii ~ 
10 aiarge "wc observc a systematic slowdown. 

The change of shape of the curves in Fig. 6(a) suggest 
that confinement induces caging dynamics. This is the 
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inhibited motion of a particle due to its "cage" of neigh- 
boring particles [47-51]. At the earliest times {At < 1 s, 
not shown), particle motion is diffusive as particles have 
not moved far enough to encounter the cage formed by 
the neighboring particles [52]. As the particle displace- 
ment becomes larger, its motion is impeded by its neigh- 
bors which form the cage, resulting in a greatly decreased 
slope of (Ax^) for At < 100 s. For smaller values of 
H, the decreasing height of {Ax'^) in this range suggests 
that the cage size decreases in more confined samples. 
This is likely due to the concentration of particles into 
the layers (Fig. 5), which crowds them within the lay- 
ers and reduces their cage sizes. Returning to Fig. 6(a), 
the upturn at larger At for (Ax^) is the result of cage 
rearrangements [7, 20, 48, 50]. The neighbors rearrange 
and this lets the caged particle move to a new position. 
The motion of particles at longer lag times is diffusive 
due to the uncorrelated cage rearrangements [50]; this is 
not quite seen in our data sets here as the time scales for 
this diffusive motion is longer than our observation times. 
The results shown are for one volume fraction; our prior 
work showed that for larger (p, the onset length scale for 
the confinement-induced slow-down increases [17]. 

To contrast the mobility in the parallel and perpendic- 
ular directions, in Fig. 6(c) we plot {x^) and (z^) sepa- 
rately for a selection of three thicknesses. Not surpris- 
ingly, the mobility is less in the z direction (perpendicu- 
lar to the wall) . Furthermore, the upturn of the MSD at 
large At is barely beginning for the z data. The contrast 
between the x and z motion suggests that cage rearrange- 
ments may favor motions parallel to the walls. 

The MSD curves show an overall slowing down due 
to confinement, but obscure the influence of the den- 
sity layers on the motion. Figure 7(a) shows the number 
density for one data set. In panels (b) and (c) we plot 
the components of the MSD, for fixed values of At, that 
are perpendicular and parallel to the walls. The dips 
in (z^) [Fig. 7(c)] coincide with the layers in panel (a) 
and imply that particles within layers are in a preferred 
structural configuration and are less likely to move else- 
where [6, 11, 17, 43, 46]. Meanwhile, the parallel compo- 
nent of motion shows no variation with z, even for long 
time scales. Our observations differ from one prior ex- 
periment by Eral et al. [11]. They found a decreased 
parallel mobility near the walls but did not measure per- 
pendicular mobility. One difference is that they stud- 
ied a single-component sample with a polydispersity of 
8%, whereas we study a binary sample. Another dif- 
ference is that their experiment had a spatial gradient 
in volume fraction due to non-density matched particles 
(they have a density difference between solvent and par- 
ticles of Ap « 800 kg/m'^, much larger than our value 
Apss 0.3 kg/m3). 

Intriguingly, our results shown in Figs. 5,7 look strik- 
ingly similar to recent experiments by Wonder, Lin, 
and Rice [41]. They studied a monodisperse quasi-two- 
dimensional colloidal system, where particles were lim- 
ited to one layer in z, and further constricted in y anal- 
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FIG. 7: (Color online) (a) Particle number-density nsmaii(2) 
as a function of distance z across the sample cell. Additional 
particles are permanently stuck to the walls of the cell (not 
shown) which have centers located a.t z = 0.00 fim and z = 
H = 8.06 ^m, indicated by the vertical dashed lines. These 
data correspond to the H = 8.06 fim data in Fig. 6. (b) Mean 
square displacement parallel to the walls (|[{Ax^) + (Aj/^)]) 
and (c) perpendicular to the walls ((Az^)) as a function of the 
particles' initial positions z. The displacements are calculated 
using a range of At, as labeled. The dotted lines indicate the 
position of the number density maximum of each layer in (a) 
while the dashed lines correspond to the approximate position 
of the centers of the particles stuck to the glass walls. 



ogous to our confinement in z. They found that their 
experimental short-time diffusion coefficients had a sim- 
ilar qualitative behavior to what is shown in Fig. 7(b,c) 
[41]. They did not study long-time diffusion coefficients. 

Our observed reduced particle mobility perpendicular 
to the walls is similar to the observations of Dullens and 
Kegel, who studied the first layer of colloidal particles at 
a smooth glass surface [28, 53]. In their work, quasi-two- 
dimensional (q-2D) layers of particles formed along the 
surface of a glass slide in a bulk polydisperse colloidal 
suspension, just as we observe. Their wall-based parti- 
cles only intermittently exchanged with the bulk parti- 
cles [53]. In their q-2D wall layer, particles exhibited 
two-dimensional behavior that was fundamentally dis- 
tinct from the dynamics of the particles further from the 
wall. However, a primary reason for this was that the 



6 



T 



0.10 
0.08 
0.06 
0.04 
0.02 
0.00 
5 
4 
3 
2 
1 

5 
4 
3 
2 
1 





1 1 


1 


1 ' 1 ' 1 ' 1 ' 1 ' 1 ' 1 






(c) 1- 


" 100 






- 50 






-25 






" 12 
~2s 











-2 2 4 6 



10 12 14 



FIG. 8: (Color online) (a) Particle number-density nsmaii(z) 
as a function of distance z across the sample cell. Additional 
particles are permanently stuck to the walls of the cell (not 
shown) which have centers located at z — 0.00 and z = 
H = 16.0 nm, indicated by the vertical dashed lines. These 
data correspond to the H > 15.8 data in Fig. 6. (b) Mean 
square displacement parallel to the walls {^[{Ax'^) + (Aj/^)]) 
and (c) perpendicular to the walls ((Az^)) as a function of the 
particles' initial positions z. The displacements are calculated 
using a range of At, as labeled. The dotted lines indicate the 
position of the number density maximum of each layer in (a) 
while the dashed lines correspond to the approximate position 
of the centers of the particles stuck to the glass walls. 



particles were fairly monodisperse, and thus could form 
monodisperse 2D phases [28, 53]. In our experiments, 
the particle layers near the wall become more pronounced 
with decreasing H, suggesting that these layers become 
more q-2D. However, our samples are binary and the two 
particle sizes remain well-mixed even at the walls (con- 
firmed by Die microscopy) . While the q-2D nature of our 
layered particles may partially explain their slow motion. 
Fig. 7 shows that slowing is not restricted to these layers 
alone. Note that the hydrodynamic interaction of parti- 
cles with nearby walls diminishes as the volume fraction 
is increased [42]. 

One explanation for the slower dynamics might be that 
the volume fraction is larger in confinement. We first con- 
sider an observation from our experiment: the pair corre- 
lation function g{r) changes slightly upon confinement, as 



shown in Fig. 9. This function indicates the likelihood of 
finding a particle a distance r away from a reference parti- 
cle at r = 0, and so the first peak position indicates a typ- 
ical spacing between nearest neighbor particles. For ideal 
hard spheres this first peak position is always at contact 
(?'max — Sagmaii). Our particlcs are slightly charged, so 
the first peak shifts to larger values. The peak is addi- 
tionally rounded by our finite resolution and the particle 
polydispersity [30]. Given the particle charges, an ap- 
proximate expectation is that (j) ~ ''max- The inset to 
Fig. 9 shows that confinement causes Tmax to shift to 
lower values, which would correspond in an increase of 
from 0.42 to 0.49. One explanation for this is that, given 
the layering of particles, the local volume fraction within 
a layer is higher than 0.42, and g{r) is reflecting this local 
volume fraction [which would be more heavily weighted 
in the average used to calculate g{r)]. Another possi- 
ble explanation is that the sample chamber is effectively 
thinner than we believe, due to the interactions between 
the particles and the walls. As noted above, the particle 
concentration is diminished within 0.4 fim of the walls. 
At the thinnest regions we study, H ^ 6 /im; if the true 
value is ~ 5.2 fim, this would correspond to an increase 
of from 0.42 to 0.42 x (6/5.2) w 0.48, consistent with 
the estimate from the g{r) data. However, an effective 
volume fraction increase from 0.42 to 0.49 seems unable 
to explain all of the dramatic slowing of the dynamics 
seen in Fig. 6(a). Consider the data at At = 100 s: 
(Aa;^) drops by a factor of ~ 40 going from the bulk to 
H — 6.25 fim. A study of an unconfined binary suspen- 
sion similar to ours found a drop of ~ 3.7 for a change of 
(p from 0.42 to 0.49 [30]. Thus we are left with a factor of 
ten in additional slowing which is not due to a possible 
volume fraction change. This agrees with the conclusions 
of Eral et al. [11]. 



C. Defining cooperatively rearranging regions 

The features of our (Ax^) curves resemble those of 
bulk supercooled colloidal liquids, where cage rearrange- 
ments play a significant role in the material's underlying 
dynamics. The process of cage rearrangements leads to 
a liquid's overall structural relaxation [54, 55]. Adam 
and Gibbs were the first to hypothesize the existence of 
"cooperatively rearranging regions" (CRRs) as a super- 
cooled liquid's means of increasing its configurational en- 
tropy [1]. Prior simulations [2, 56, 57] and experiments 
[20, 58, 59] found cooperatively moving regions, defined 
as groups of neighboring molecules or particles that col- 
lectively rearrange their positions. The connection be- 
tween these observations and the CRRs of Adam and 
Gibbs is perhaps problematic [60] and there are other 
theories that also consider CRRs [3, 57, 61]. Independent 
of the theoretical status of CRRs, it is certainly intrigu- 
ing that spatially heterogeneous dynamics have been seen 
in a wide range of glass-forming systems [3] . We wish to 
see how the character of spatially heterogeneous motions 
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FIG. 9: (Color online) The pair correlation function g{r) for 
a range of H. Darker curves correspond with thinner samples. 
The curves are from samples with thickness H equal to 6.25, 
6.63, 8.09, 12.6, 15.8, and 18.9 fim, along with one curve 
for the sample's bulk (the lightest color curve). The inset 
shows the position of the first peak as a function of H. The 
red horizontal dashed line indicates rmax for the bulk sample, 
while the diagonal dotted black line is a guide to the eye. 



changes upon confinement. 

The precise definition of a cooperatively rearranging 
region is open to interpretation. Our definition is de- 
scribed below, and is comprised of three key elements: 
(1) the time scale used to determine displacements, (2) 
the threshold for considering a displacement to be a "re- 
arrangement," and (3) the definition of which particles 
are adjacent, such that their motion is "cooperative." 

We first define the time scale of interest. Prior work 
found that a good choice is based on the shape of 
the probability distribution of displacements. Rearrang- 
ing particles have displacements which are larger than 
normal, and thus lie in the tails of the distribution 
[20, 58, 59, 62, 63]. The size of the distribution tails 
is quantified by the non-Gaussian parameter a2, 



a2(Ai) 



3(Aa;2 



1, 



from Ref. [64] . The maximum of a2 defines the cage rear- 
rangement time scale At* . We plot a2{At) in Fig. 6(b): 
both the maximum value of a2 and the time scale A<* 
increase with decreasing H, similar to prior observa- 
tions on a monodisperse sample [19]. For data from 
H > 15.8 /im, the levels of noise at low values of At 
manifest as a false increase of 012, so we ignore this peak. 
For At > 10 s there are secondary local maxima of 02 
that we consider to be a better determinant of At* [see 
the inset plot of Fig. 6(b)]. We plot At* versus H in 
Fig. 10(a), which decays roughly exponentially with H 
until iJ 20 /xm, at which it reaches the bulk value. 
Simply put, as H decreases the displacement distribu- 
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FIG. 10: (Color online) (a) Values of At* that maximize a2 
for a range of thicknesses H. The dotted line is a guide to 
the eye. The horizontal dashed line indicates the value of At* 
for the bulk sample, (b) Plot of the mean number of particles 
within cooperatively rearranging regions, as a function of H. 
Only Nc > 3 are considered to avoid trivial rearrangements 
that consist of 1 or 2 particles. The upper dashed line indi- 
cates (Nc) for the sample's bulk, (c) Plot of the average ex- 
tent of cooperatively rearranging regions parallel (crosses and 
circles correspond to x and y, respectively) and perpendicular 
(triangles, z) to the sample's walls. The difference between 
the X and y data is an indication of the amount of uncertainty 
in our data. The upper and lower horizontal dashed lines in- 
dicate the mean horizontal and perpendicular extent of data 
from the sample's bulk, respectively. In principle these should 
be the same (the behavior should be isotropic in the bulk); in 
practice the lines may differ due to finite data or anisotropy 
in the imaging volume (50 fim in x and y but only 20 fim 
in z). (d) Value of the mean square displacement (Ar^) at 
the time scale At=100 seconds. The horizontal dashed line 
indicates the value of the bulk sample. 



tions become less Gaussian-like, and the time scale At* 
for which the distributions are most extreme grows. 

To define the length scale which separates a "rearrang- 
ing" displacement from a "caged" displacement, we use 
a mobility threshold Ar*. Both experiments [20, 65] and 
simulations [66] have used a displacement threshold to 
define mobility such that over time, some percentage of 
the particles have displacements |Af^ > Ar* [2, 66], al- 
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though at any given time the fraction may not be exactly 
this percentage. Thresholds of the top 5**^ percentile 
[20, 66], 8*^ percentile [50], lO*'' percentile [65], and 20*'' 
percentile [67] have all been used to define Ar*. From ex- 
amining distributions of Ax and Ay for our data for each 
H, we find that the slowest 90% of the displacements arc 
well described by a Gaussian distribution, whereas the 
top 10% are more probable than a Gaussian distribution 
would predict. Thus, we define our mobility threshold as 
the top 10% of the most mobile particles. Displacements 
in the z-direction, however, vary significantly with i?, 
making their inclusion in the calculation of our threshold 
impractical. As is the case with prior studies [2, 20], our 
choice of Ar* is somewhat arbitrary and our results are 
robust to some variation of Ar*. 

To complete our identification of CRRs we must iden- 
tify which highly mobile particles are simultaneously 
nearest neighbors. Similar to other work, we define neigh- 
bors as those particles whose separation is less than a 
cutoff distance set by the first minimum of the pair cor- 
relation function g{r) [50]. Our distributions of g{r) do 
not vary substantially with H, as shown in Fig. 9. We 
use the average position of the first minimum (3.87 /im) 
to define particles which are nearest neighbors. 

One problem we face is the selective visibility of the 
colloidal particles. As discussed earlier, only the smaller 
particles of our binary suspension are fluorescently la- 
beled meaning that the larger species of particles are not 
visible to our confocal microscope. Despite this limita- 
tion we can still draw some reasonable conclusions. For 
example, in a study of the aging of a binary colloidal 
suspension similar to the one studied here. Lynch et al. 
showed that the cooperative dynamics of one species were 
similar to that of the other [67]. Mobile particles of one 
species were usually near mobile particles of the other 
species. Therefore, it is reasonable to draw some conclu- 
sions about cooperative motion from the small particles 
alone. One other related limitation is that small rear- 
ranging particles may not be nearest neighbors, but may 
be part of the same CRR, connected by unseen large par- 
ticles. This is beyond our ability to determine, although 
it may simply limit the apparent sizes of CRRs without 
otherwise changing their character. 



D. Shapes of Cooperatively Rearranging Regions 

We first visualize these CRRs to develop a qualitative 
understanding of their nature. Figure 11 depicts clusters 
of the top 10'^ percentile of the most mobile particles in a 
sample confined within a plate-spacing oi H = 15.8 /im 
and 6.63 /im (panels (a,b) and (c,d) respectively). For 
clarity, bonds have been drawn between particles that 
are nearest neighbors, i.e. within a cluster. For both 
thicknesses, groups of mobile particles can be seen. The 
size of these mobile clusters in the unconfined sample is 
small, as expected for this low volume fraction (0 = 0.42) 
[20]. Despite their small size, these mobile clusters are 



the primary means for particle rearrangments in the sam- 
ple. The sample can be considered as composed of the 
slowest 90% particles which are caged at a given mo- 
ment, and the rearranging fastest 10%. If the nature of 
the fastest 10% changes in confinement - for example, 
if those rearrangements occur less frequently - then the 
overall sample diffusivity will decrease. 

Confinement induces slower dynamics, and in the bulk 
slower dynamics are associated with larger CRRs [20, 59]. 
Perhaps confinement induces a similar larger size of 
CRRs [19]; but at first glance, comparing Fig. 11 pan- 
els (a) and (c) might suggest that the cluster sizes are 
smaller upon confinement. However, recall that the par- 
ticles shown are the most mobile 10%; the thinner sample 
has fewer particles in the imaged volume, and thus 10% 
of this smaller number results in fewer mobile particles 
to show without necessarily implying that the CRRs are 
smaller. To quantify the size of CRRs we calculate the 
mean number of particles in a CRR Nc as a function of 
ff, plotted in Fig. 10(b). Figure 10(b) shows that CRRs 
tend to involve roughly the same number of particles, re- 
gardless of thickness. The mean CRR size is between 3 
and 4 particles, but this is only slightly larger than the 
minimum size of 3 particles. The small size may be be- 
cause the bulk sample, with (p = 0.42, is liquid-like and 
only has small CRRs [20]. Alternatively, as noted above, 
we cannot see the large particles which are almost cer- 
tainly part of CRRs [67]. With the data of Fig. 10(b), we 
cannot say clearly if the CRRs are larger or smaller upon 
confinement. There is a very slight downward trend in 
{Nc) with decreasing H, but this could be due to poor 
statistics. It is possible that the influence of confinement 
on the size of CRRs would be clearer in a sample with a 
larger value of although such samples are very difficult 
to load into our thin sample chambers (as has been noted 
by others [68, 69]). Likely some of the difficulty in loading 
the samples is due to their increasing glassiness in con- 
fined spaces. Results from another confocal microscopy 
experiment on a monodisperse sample suggested that the 
length scale for CRRs grows upon confinement [19]. The 
difference from our results may be due to our use of a 
binary sample. 

An alternate way to quantify the size of a CRR is 
through its spatial extent. We define the spatial extent 
of the CRRs as ^extent — max(a;j) — min(xi), where i 
ranges over all particles within a given cluster of mobile 
particles. Similar definitions apply for the y and z di- 
rections. We plot the mean CRR extent in the x, y and 
z directions separately in Fig. 10(c). We find that the 
CRRs maintain a constant size in the direction parallel 
to the walls. However, the amount of distance that the 
CRRs extend in the direction perpendicular to the walls 
is significantly smaller than i?, and decreases as H de- 
creases. In the z direction, then, clusters are smaller, 
perhaps trivially because CRRs have to fit into a thinner 
sample chamber. In the x and y direction, clusters main- 
tain a constant size with confinement. Comparing this 
result with the (TVc) data of Fig. 10(b) suggests that the 
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FIG. 11: (Color online) Three-dimensional renderings of the top lO*'' percentile of the most mobile particles in two different 
sample thicknesses. The gray (green online) bonds between particles are drawn only to indicate nearest neighbors and do not 
imply permanent attachment between the particles. Only groups with Nc > 2 particles are drawn. The colors correspond 
to the magnitude of displacements in the z-direction, normal to the confining boundaries. The experiment's field of view, and 
the effective position of the confining boundaries, are indicated by the light gray bounding boxes. The sample on the left (a, 
b) has a thickness of H — 15.8 fim and the one on the right (c, d) has H — 6.63 fj,m. The top row of images (a, c) view 
the sample normal to the confining boundaries, while the bottom row (b, d) provide a parallel view. Black and white indicate 
displacements of at least asmaii over a At = 23 s and 250 s for the thicker and thinner sample respectively. 



CRRs are becoming more compact in z with the same 
number of particles. This suggests that perhaps they are 
fractal in the bulk with a fractal dimension larger than 2 
(as seen previously in Ref . [20] ) and become more planar 
upon confinement (fractal dimension approaching 2). 

The onset of flatter or more planar CRRs coincides 
with the sample's overall slowing. In Fig. 10(d) we plot 
the MSD values from Fig. 6(a) for At — 100 s against the 
corresponding range of H. We observe that the MSD val- 
ues of Fig. 10(d) begin to deviate from those of the bulk, 
indicated by the horizontal dashed line in (d) , at approx- 
imately the same H that the z-extents of the CRRs first 
begin to flatten relative to the z-extent from the bulk, 
the horizontal dashed line in (c). This is the strongest 
evidence linking the changing CRRs to the slowing dy- 
namics. The overall concept is that confinement modi- 
fies the structure from that of the bulk, and this changed 
structure leads to slower dynamics [45]. 

We can more carefully quantify the shapes of the CRRs 



by considering the probability distributions of the extents 
in the three directions. These distributions are shown in 
Fig. 12 for a bulk sample (panel a) and a confined sam- 
ple (panel b). In the unconfined sample the probabilities 
of the extent in the x, y, and z directions are approxi- 
mately the same, as should be expected; these CRRs are 
spatially isotropic. Differences in the z are most likely 
due to minor particle position errors which are larger 
in z, as discussed in Sec. II. In contrast with Fig. 12(a), 
Fig. 12(b) shows that the extents of CRRs in confinement 
have a very different probability distribution. The extent 
in z is nearly zero for a majority of CRRs [red (dark gray) 
curve in Fig. 12(b)]; these are planar CRRs and are over- 
whelmingly more probable than in the unconfined case. 
A small subset of confined CRRs do extend into the z- 
direction by one to two particle diameters. The clusters 
of rearranging particles along the walls in Fig. ll(c, d) 
seem to be the most planar in shape. Thus we are led to 
conclude that the CRRs in the confined cases are quali- 
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FIG. 12: (Color online) Extent of mobile groups of parti- 
cles. Dark and light gray (red and orange online) indicate the 
perpendicular and parallel extents respectively. CRRs in the 
bulk (a) of the sample are more isotropic in shape than when 
confined (b) to H = 6.63 /im. Note the difference in vertical 
scales used by the two plots. Only CRRs consisting of at least 
3 particles are considered in the data shown. 



tatively different than those of the unconfined sample. 



E. Details of Rearrangements 

We next investigate the behavior of particles within 
CRRs. In Fig. 11 the particles are colored in corre- 
spondence with their amount of perpendicular motion, as 
shown in the key. In the confined situation mobile parti- 
cles displace horizontally more frequently than otherwise, 
as suggested by the greater number of orange (medium 
gray) particles in Fig. ll(c, d). This makes sense: a re- 
arrangement consisting of particles within a single layer 
does not require the particles to move vertically for the 
rearrangement to occur. Occasionally we do see particles 
which jump between layers or even swap between layers; 
one example is near the bottom right corner of Fig. 11(c). 

To compare the amount of parallel versus perpendicu- 
lar displacements, we calculate the directions of motion 
for all particles and then repeat the comparison for dif- 
ferent confinement thicknesses. Using a spherical coor- 
dinate system we determine the polar angle of a given 
particle displacement. The polar angle 6 spans a range 
from 0° to 180°, which correspond to motion toward or 
away from the nearest sample chamber wall, respectively. 
That is, we exploit the symmetry between the two walls. 
We first compute the polar angle 9 relative to the +z axis, 
and then use 180° — 9 for the data in the lower half of the 
sample chamber. Comparing the data separately for the 
top and bottom half, we find no difference in the results. 
For isotropic motion, the distribution of 9 is proportional 
to sin 9, so we divide our measured histograms by sin 9 



FIG. 13: (Color online) Polar plots of the probability 
distributions of the directions of particle displacements for 
H = 6.63 and H — 15.8 as indicated. (a,b) Data for 
the two thicknesses, considering only particles away from the 
walls. (c,d) Data for particles in the layers immediately adja- 
cent to a wall. The displacements from one wall are reversed, 
so that 180° always means motion away from the nearest wall. 
In all panels, the light blue curve (light gray) is the distribu- 
tion for the most mobile 10% of the particles, while the dark 
blue curve (dark gray) is the distribution for all particles. Dis- 
placements are measured over At = 250 s and At = 23 s for 
the H = 6.63 /xm and H = 15.8 jum data, respectively. 



to remove this dependence. The distributions are plotted 
in polar coordinates, shown in Fig. 13, for thicknesses of 
H = 6.63 /im and H = 15.8 /im. The dark curves are 
for all particles, and the light curves are for the top lO*'' 
percentile of displacements, providing insight into the di- 
rections that tend to permit higher mobility. The top 
panels show the motion of the particles in the interior of 
the sample, and the bottom panels show the motion of 
the particles immediately adjacent to the walls. 

In both the 15.8 /im and 6.63 fim samples, the particles 
in the outer layers along the walls tend to move parallel 
to them {9 90°) rather than perpendicularly [Fig. 13(c, 
d)]. The effect is even more pronounced for the fastest 
particles, whose distribution suggests that fast particles 
move almost exclusively along the walls. This agrees with 
our observations from Fig. ll(b, d), where the particles 
layered along the walls are almost all orange (medium 
gray) , indicating they are moving primarily horizontally. 
The distributions in Fig. 13(c, d) do show some data at 
9 — 180°, indicating that some particles move away from 
the walls, and less data at 6* = 0°, indicating that some 
particles make slight motions toward the walls. 

The situation changes markedly for the inner layers 
[Fig. 13(a, b)]. Considering only the full distribution 
of all particles we see that the displacements are more 
isotropic, although there is still a slight bias in the 
9 90° direction. The distribution of directions for the 
most mobile interior particles is similar. There are bumps 
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in these distributions near 6* = 0° and 180°, which sug- 
gests that particles that move in z have a slight increased 
probability to make large motions in z, hopping between 
layers. 

Overall, the particle dynamics in the thicker region are 
far more isotropic than the ones from the confined region 
[compare Fig. 13(a, c) with (b, d)]. In the H = 15.8 /im 
case, there are appreciable signs of anisotropic behavior 
only along the walls. 

IV. CONCLUSION 

The smooth quasi-parallel walls confining our sample 
induce the formation of density layers within the col- 
loidal sample's volume. The most dense layers form along 
the sample chamber's glass surfaces, as shown in Fig. 5 
and also observed in other experimental work that used 
single-component colloidal samples [11, 28]. The struc- 
tural inhomogeneities induced by the density layers result 
in corresponding inhomogeneities in the system's dynam- 
ics, as described by the plots in Fig. 7. Particles move 
most easily within their layer, but this is still slower than 
they would move in unconfined samples. The layered par- 
ticles cooperatively rearrange within the layer but rarely 
with adjacent layers; the cooperative rearrangements oc- 
cur in more planar-shaped groups of particles. Given that 
even in unconfined samples, particles need to move co- 
operatively if they wish to have large displacements, the 



change in the character of the cooperatively rearranging 
regions seems to explain the slowing dynamics. In short, 
the thickness at which we begin to observe the slowing in 
the sample's average dynamics corresponds with the con- 
finement length scale at which cooperatively rearranging 
regions begin to become planar in shape [Figs. 10(c) and 
(d)]. Our prior work suggests that the observed increase 
in rearrangement time scales and the thickness at which 
these regions begin to flatten will both grow with higher 
volume fractions [17]. 

It is likely if the walls were roughened, the results 
might change. Simulations [6, 7, 10] and experiments 
[11, 36] showed that behavior is often glassier with rough 
walls. With rough walls, layering is greatly diminished or 
prevented entirely, or perhaps becomes more subtle. For 
example, particles might form a corrugated layer wrap- 
ping around the local wall texture. This could then lead 
to other shapes for the cooperatively rearranging regions; 
the main point being that structure that departs from the 
bulk results in slower dynamics [45]. 
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